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We update our exploration of the MSSM parameter space at the weak scale 
where new accelerator and cosmological constraints are respected. The depen- 
dence of WIMP-nucleon cross sections on parameters of the MSSM, uncertainties 
of the nucleon structure and other theoretical assumptions like universality and co- 
annihilation are considered. In particular, we find that the coannihilation does not 
^ . have a significant effect in our analysis in certain regions which are allowed even 

. with coannihilation. The new cosmological constraint on the relic neutralino density 

1^ ( , used in the form 0.1 < ^Ix^o < 0-3 does also not significantly affect the regions of 

allowed neutralino-nucleon cross sections. We notice that for nuclear targets with 
spin the spin-dependent interaction may determine the lower bound for the direct 
detection rate when the cross section of the scalar interaction drops below about 
10-12 pb. 
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I. INTRODUCTION 



' It is well known (see, for example, that the Minimal Supersymmetric Standard Model 

, (MSSM) being the most promising extension of the Standard Model offers a solution of the 

r*^ ' hierarchy problem, possesses gauge coupling unification and naturally proposes a Dark Matter 

C^i candidate — the lightest supersymmetric particle (LSP). In the framework of the low-energy su- 

persymmetry (SUSY), when SUSY breaking masses lie below a few TeV, sparticles will be copiously 
produced (and detected) at future coUiders like the Large Hadron Collider (LHC) at CERN. On 
. the other hand, there are several on-going and future projects searching for the LSP as a Dark 

Matter particle. One of them even claims a positive signal ||], although the situation remains 
rather contradictory . The present experimental upper limit on the spin- independent part of the 
^ ■ elastic scattering of the LSP on a nucleon is around 10"^ pb for 50 GeV < < 100 GeV. In any 

^ I case, it seems very plausible that both, SUSY collider signals and LSP Dark Matter will be found 

in future. Such dark matter searches offer interesting prospects for beating accelerators in the 
discovery of SUSY, particularly during the coming years before the LHC enters operation Q. In 
this situation naturally arised the question of how small the event rate of the LSP direct detection 
can be, provided the LSP is a cold dark matter particle. Searching for the answer different SUSY 
models were considered (see for example, [||-||). 

Recently exciting evidence for a flat and accelerating universe was obtained [|l^,0 . The position 
of the first acoustic peak of the angular power spectrum strongly suggests a flat universe with 
density parameter fJo = 1 while the shape of the peak is consistent with the density perturbations 
predicted by models of inflation. The density parameter f2o = Po/Pc is the ratio of the current 
mass density po to the critical density pc — 1.88 • 10~^^ hg g-cm"^, with fen b eing the dimensionless 
Hubble parameter. Data support the straight line flo = flM + = 1 Il2|-|l5|], where f2M is the 
matter density in the universe and Q\ is the contribution of the non-zero cosmological constant 
(the energy density of the vacuum). At the same time one determines fi^/ = 0.4 ± 0.1, which 
implies fl\ = 0.85 ± 0.2, a value that has been supported from high-redshifted Supernova data 
||lq| . Since the baryonic matter density is small, fls — 0.05 ± 0.005, the values for matter density 
flM give a Cold Dark Matter (CDM) density JIcdm — 0.35 ±0.1, which combined with recent 
measurements of the Hubble parameter ho — 0.65 ± 0.05, result in smaller CDM relic densities 
ncBM feo' 0.15 ±0.07 0-|l|. 
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Previously we have restricted our analyses to the cosmological constraint for the relic density 
of the LSP in the range 0.025 < ^x^o ^ ^ lOHl^]' where the neutralino relic density parameter 
^ Pxl ^'^d Px is the relic neutralino mass density. In this paper adopting the above-mentioned 
new cosmological data and going to compare our estimation with literature y,^,^,^ we constrain 
the allowed region for the neutralino relic density in the form 0.1 < ^x^^ < ^.3. It is possible that 
there is more than one component in the cold dark matter, so that < ^^cdm, and therefore 
< 0.1. Although, in general, a decrease of Vl is associated with larger elastic scattering cross 
sections, the detection rate also must be reduced because of the corresponding reduction in the 
density of LSPs in the Galactic halo. Here we neglect this possibility and assume that all the cold 
dark matter is composed of LSPs, so that fi^ > 0.1 

There are two main approaches to evaluate the dark-matter-neutralino nucleon cross section and 
the expected event rate in a detector. The basis of the first approach is the minimal supergravity 
(mSUGRA) model | p3[ . This model assumes the minimal supersymmetric standard model to be 
valid at all energy scales from Aiwcak up to A/qut — 2 x 10^^ GeV. The mSUGRA model arises as 
the low-energy limit of a supergravity theory, where supersymmetry is broken in the hidden sector 
of the model at an energy scale of M ^ 10^° GeV. Supersymmetry breaking is communicated to the 
observable sector via gravitational interactions, leading to soft SUSY breaking mass terms of the 
order of the electroweak scale. At the GUT scale this leads to a common mass for all scalars mo and 
a common trilinear coupling Aq. Motivated by the apparent unification of gauge coupling constants, 
it is also assumed that all gaugino masses are unified to mi/2 at AfcuT. The weak scale sparticle 
spectrum is derived from renormalization group running of the SUSY soft breaking parameters. 
Requiring radiative electroweak symmetry breaking allows the determination of the superpotential 
Higgsino mass squared /i^, and allows the expression of the soft SUSY breaking bilinear term B 



clings are derived in terms of the minimal 
-E5|. From a practical point of view this 



in terms of tan/3. Thus, all sparticle masses and con 
parameter set toq, ^1/2; ^o, tan/3, and sign(/i) |^ 
approach as much as possible relies on theoretical arguments like unification, naturalness etc, 
aiming to maximally reduce the set of free parameters and obtain maximally restricted predictions. 
In this approach Q the quantum stability of the gauge hierarchy suggests that sparticles weigh 
less than about 1 TeV ||2^, which is also the range favored for a cold dark matter particle, and 
there are indeed generic domains of the MSSM parameter space in which the relic LSP density falls 
within the range favored by astrophysics and cosmology. The unsuccessful laboratory searches for 
sparticles impose non-trivial constraints on the MSSM parameter space, suggesting that the LSP 
neutralino is mainly a C/(l) gaugino (Bino) In the MSSM the lightest neutralino x = Xi is a 
mixture of four superpartners of gauge and Higgs bosons (Bino, Wino and two Higgsinos): 

X = Nix&' + Ni2W^' + Ni:iH^I + NuHl (1) 

It is commonly accepted that x is mostly gaugino-like if P = N^i + /Vfj > 0.9 and Higgsino-likc if 
P < 0.1, or mixed otherwise. 

It is due to the Bino-likeness of the relic neutralinos that the calculated neutralino-nucleon cross 
sections appeared very small and one usually arrived at the conclusion that it was hardly possible 
to reach the mSUGRA space by means of direct and indirect searches for dark matter particles 
P,p|,p|, pj| , p4[[2^ pO| . The other side of this conclusion is also well known: when decreases, the 
Higgsino components /V13 and 7Vi4 of x increase (P decreases) and as a result the spin-independent 
cross section increases. So Higgsino-like (and mixed) neutralino on the other hand increases the 
prospect for its detection as a dark matter particle |7[ |2^ , |3l| , |3^ . Therefore it seems a crucial 
question here, to what extent the neutralino is mostly gaugino-like, Higgsino-like, or mixed. 

A way to look for any possibility of higher cross sections and higher expected rates of detection 
is to investigate alternate models. The basis of it is a departure from the stringent mSUGRA by 
means of a relaxation of some unification and other theoretical assumptions aiming to obtain as 
general predictions for the expected detection rate as possible. A remarkable shift from mSUGRA 
to more relaxed models was made by . It mostly included relaxation of unification of soft 

scalar mass parameters (so called nonuniversal soft symmetry breaking) as well as gaugino mass 
non- universality. The large tan/3 regime was also considered as a source of higher cross sections. 
Indeed in canonical mSUGRA it was pointed out ||3^-|35|| that the large tan j3 regime allows regions 
where u^p ~ 10~^pb. Besides, with non-universal soft scalar masses, it was also found that 
cTxp ~ 10~^ pb for small values of tan/3. In particular, this was obtained for tan /3 > 25 (tan/3 > 4) 
working with universal (non- universal) soft terms in [ p5| . These analyses were performed assuming 
(non-)universality of the soft breaking terms at the unification scale, Mqtit ~ 10^^ GeV, which 
can be obtained within superstring theories Pq,p3l ^-nd heterotic M-theory |3q,pq| . 
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Completely new possibilities have also been discussed. For example, it is found that in super- 
symmetry models muli-TeV scalar masses can exist consistent with naturalness on a certain branch 
of the radiative breaking of the electroweak symmetry . A similar phenomenon appears in the 
so-called focus point supersymmetry models where one also avoids naturalness constraints with 



multi-TeV scalars |31 4^,^, and in models with moving intermediate unification scale [ p2[ . 

It was noticed that the assumptions concerning universality of the scalar masses mi(A/GUT) = 
mo, and the trilincar scalar couplings v4''"''^(Mgut) = ^o, are not very solid, at least from phe- 
nomenological point of view, since, universality might occur at a scale higher than Mqut 
10^^ GeV Q, or according to string models at a scale M/ smaller than Mgut ^ GeV 



It was realized that the string scale may be anywhere between the weak scale and the Planck 



scale. For instance D-brane configurations allow these possibilities in type I strings |44-47|. Sim- 



ilar results can also be obtained in type II strings [j48| and weakly and strongly coupled heterotic 



strings 1 49 5^] . Moreover the M/ might be anywhere between the weak scale and the Planck scale 
, with significant consequences for the size of the ncutralino-nucleon cross section. The case of 
non-universal gaugino masses was analyzed in |^,^^J3^ and with respect to direct detection of 
the superlight dark matter neutralinos in jiot . Schemes with CP-violating phases one can find in 

il- 

Therefore due to the large uncertainties involved in the choice of the scale Mj and going to obtain 
as much as general predictions it appeared more convenient to work within a phenomenological 
SUSY model whose parameters are defined directly at the electroweak scale as for example in 
§,|,||-|6| and which is denoted as an effective scheme of MSSM (effMSSM) in ||. 

Obviously, this way much larger expected event rates were obtained and optimistic conclusions 
concerning the possibility to constrain significantly the SUSY parameter space with dark matter 



experiments were drawn |8,P4t51,R2 53 



In our previous calculations in effMSSM P,p_8 19 5q] we have adopted an (effective) scheme (with 



non-universal scalar masses and with non-universal gaugino soft masses) which has supplied us with 
large relatively direct detection rates of dark matter neutralinos, practically independently of what 
is the neutralino composition. In most of the MSSM parameter space we, in agreement with 
others, have obtained at the detectable level mostly gaugino-like neutralinos, but always existed 
small Higgsino admixtures (at a level less than 1-5%) which managed to produce large enough 
cross sections and rates. 

In 1994 we claimed that nuclear spin is not important for detection of dark matter particles, 
provided the detection sensitivity does not exceed 0.01 events/day/kg, which was considered that 
time as unreachable |5|]. Now the situation has changed and we would like to notice that for 
targets with spin-non-zero nuclei it might be the spin-dependent interaction that determines the 
lower bound for the direct detection rate when the cross section of the scalar interaction, which is 
usually assumed to be the dominant part, drops below lO^^^^^'^pb 

New updated parameters of the nucleon structure involved in the evaluation of the elastic neu- 
tralino nucleon scattering have become available Q and one expects that they will affect the cross 
sections of neutralino nucleon scattering. At least significant cancellations may occur for some 
values of tan /3 for scalar- and spin-dependent cross sections (at least for tan/3 < 10 j|]). 

The above considerations stimulated us to perform a recalculation of our previous analysis. 



II. APPROACH 



A dark matter event is elastic scattering of a relic neutralino from a target nucleus producing 
a nuclear recoil which can be detected by a suitable detector The differential event rate in 

respect to the recoil energy is the subject of experimental measurements. The rate depends on 
the distribution of the relic neutralinos in the solar vicinity and the cross section of neutralino- 
nucleus elastic scattering. In our analysis we use the total event rate R which is integrated over 
recoil energies and useful for searching for domains with extreme rates. We follow our papers 
jl^Jl^, where one can find all relevant formulas and astrophysical parameters. We consider only 
a simple spherically symmetric isothermal distribution and do not go into details of any possible 
uncertainties (and/or modulation effects) of the Galactic halo WIMP distribution ||60|-|6^. 

To calculate the event rate we use for the relic neutralino mass density and for the escape 
neutralino velocity commonly accepted values 0.3 GeV/cm'^ and 600 km/s, respectively. Their 
experimental variations can slightly change R but leave the dependence of R on the MSSM param- 
eters unaffected. To compare our results with other calculations and sensitivities of different dark 
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matter experiments we calculate also the total cross section for relic neutralino elastic scattering 
on the nucleon. The scalar (spin-independent) part of the elastic neutralino-proton(neutron) cross 
section at zero momentum transfer g = is 

agf(0)=4^[Cp,„]^ where C,,„ = 4^^^, + A/^^^^ ^ C,. 

q—u,d.s c,b,t 

The spin-dependent part of the elastic x-nucleon cross section can be written as 
CTPj^(O) =4^ 3 Up J' where ^p^„ = ^ ^^Ag^P-"), and 



The effective couplings Aq and Cq of the neutralino-quark Lagrangian 

LcS = Aq ■ X7m75X • 97^759 + CgXX -qq + O (l/m|) 
which enter the cross sections one can find in [Q. The parameters f^q'^^ and Z^^""* are defined by 

q—u^d,s 

Following Q we use the updated parameters: 

= 0.020 ± 0.004, f^Pj = 0.026 ± 0.005, f^] = 0.118 ± 0.062; (2) 
f^'J = 0.014 ± 0.003, f^'J = 0.036 ± 0.008, /^"^ = 0.118 ± 0.062. (3) 

Our estimations of the effect of the inaccuracy in the determination of fj.^ on the total event rate 
agree with those obtained before in ||5^ and in |^,^ 67 6^. For a different determination using 



an analytic analysis sec |^]. The two corridors do ovelap. The inaccuracy maximally changes the 
proton-neutralino cross section (event rate) within about one order of magnitude. The value chosen 
in this work gives probably a more pessimistic view of the cross sections. The inaccuracy of other 
parameters has a smaller effect on the cross sections. 

The factors A^-'''"-' parameterize the quark spin content of the nucleon. A global QCD analysis 
for the gi structure functions including 0{al) corrections, corresponds to the values [Q 

^(P) ^ a|;") ^ 0.78 ± 0.02, A'^f = A|;') = -0.48 ± 0.02, A^p^ ^ A^") = -0.15 ± 0.02. (4) 

We calculate ^x^o following the standard procedure on the basis of the approximate formula 
||7^,|7l|. We take into account all channels of the x ~ X annihilation. Since the neutralinos are 
mixtures of gauginos and higgsinos, the annihilation can occur both, via s-channel exchange of the 
Z° and Higgs bosons and t-channel exchange of a scalar particle. This constrains the parameter 
space |2^ , [70| |. As mentioned in the introduction we require 0.1 < flx^o < O-^; for comparison we 
also present our results when 0.025 < ^x^a < 1- 

Another stringent constraint is imposed by the branching ratio of 5 — > 37 decay, measured by the 
CLEO collaboration to be 1.0 x 10"-* < B(6 ^ 57) < 4.2 x 10"-*. In the MSSM this flavor-changing 
neutral current process receives contributions from H^-t, x^-i and g-q loops in addition to the 
standard model W-t loop. This also restricts the SUSY parameter space |Q. 

The masses of the supersymmetric particles are constrained by the results from the high energy 
colliders. This imposes relevant constraints on the parameter space of the MSSM. In |^ we used the 
following lower bounds for the SUSY particles pSl: M-+ > 65 GeV for the light chargino, M~+ > 
99 GeV for the heavy chargino, M^o ^ ^ > 45, 76, 127 GeV for non-LSP neutralinos, respectively; 
Mc, > 43 GeV for sneutrinos, Mg^ > 70 GeV for selectrons, Mq > 210 GeV for squarks, M^^ > 
85 GeV for light top-squark, M^o > 79 GeV for neutral Higgs bosons. Men > 70 GeV for the 
charged Higgs boson. On the basis of last LEP results p3 we use now new lower limits for 
charginos: M~± > 100 GcV, and neutral Higgs bosons: m„o > 100 GeV. 

As previously |^ we explore the MSSM parameter space at the weak scale relaxing completely 
constraints following from any unification assumption. Nevertheless we respect other available 
restrictions from cosmology, accelerator SUSY searches, rare FCNC b ^ decay, etc [|o|j|j2^. 
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The MSSM parameter space is determined by entries of the mass matrices of neutralinos, 
charginos, Higgs bosons, sleptons and squarks. The relevant definitions one can find in The list 
of free parameters includes: tan/3 is the ratio of neutral Higgs boson vacuum expectation values, 
jj, is the bilinear Higgs parameter of the superpotential, Mi^2 are soft gaugino masses, Ma is the 
CP-odd Higgs mass, m~, m~, m~ ('ti—, '™~) are squared squark (slepton) mass parameters for 

the 1st and 2nd generation, m~ , m~, m~ {m~ , mi.) are squared squark (slepton) mass param- 

Q3 T B L3 r 

eters for 3rd generation and At, At, Ar are soft trilinear couplings for the 3rd generation. In our 
numerical analysis the parameters of the MSSM are randomly varied in the following intervals: 

- 1 TeV < Afi < 1 TeV, -2 TeV < M2, ^i,At<2 TcV, 

1 < tan/3 < 50, 60 GeV < Ma < 1000 GeV, (5) 
10 GeV^ < ml,ml,ml^,ml^ < 10^ GeV^ 

Following |ll|,p|j67| we assume that squarks are basically degenerate. Bounds on flavor-changing 
neutral currents imply that squarks with equal gauge quantum numbers must be close in mass. 
With the possible exception of third generation squarks the assumed degeneracy therefore holds 
almost model- independently Q. Therefore for other sfermion mass parameters we used the rela- 
tions TO— — m?~ — TO— , TO— = m? , TO— = TO— = TO?, , TO— = TO? . The parameters Ah and 

U D Q E rp g Q3' Ea ^ 

At are fixed to be zero. We consider the domain of the MSSM parameter space, in which we 
perform our scans, as quite spread and natural. Any extra expansion of it like, for example, using 
— 10 TeV< M2 < 10 TeV, etc, of course, can be possible, but should be considered, contrary to 
||3^,^,0, as quite unnatural. 



III. RESULTS AND DISCUSSION 



A. Coannihilation 



The effects of coannihilation may become important when the next to the lowest supersymmetric 
particle (NLSP) has a mass which lies close to the LSP mass jTsj. The size of the effects is 
exponentially damped by the factor e"^'^ where — {rrii/m^ — 1), x = m^/kT and where to^ 
is the LSP mass. Because of this damping the coannihilation effects are typically important only 
for regions of the parameter space where the constraint < 0.1 is satisfied. Some of the possible 
candidates for NLSP are the light stau fi, en, the next to the lightest neutralino X2: and the light 
chargino Xi ■ It was found that in mSUGRA the upper limit on the neutralino mass consistent 
with the current experimental constraints on the relic density is extended from 200 GeV to 600 
GeV 1?^] when the effects of x — f coannihilation are included. 

By means of excluding points which can give non-negligible contribution to relic neutralino 
annihilation via coannihilations with other SUSY particles we simply have estimates of the influence 
of the coannihilation to our previous results. Wc used the constraint: (m^ — mi^gp) / (m-^gp) < 0.2, 
where i runs over ncxt-to-LSP neutralino, charginos, staus, stops etc. In accordance with previous 
estimates of |3|,|^,||,|l| we found that the coannihilation does not significantly change our main 
results. In fact, less than 20% of the models were denied by this coannihilation constraint, which 
in the case of 0.025 < O^/iq < 1 excludes points with simultaneously small |/i| < 500 GeV) 
and large \Mi\ {\Mi\ > 600 GeV), allowing for a substantial non-gaugino fraction of the LSP only 
in the region of relatively small If the relic abundance of neutralinos is located in the range 
0.1 < ri^/iQ < 0.3, the coannihilation constraint appears less restrictive. 



B. Cross sections 

Our calculations with the updated nucleon structure [Q for the WIMP-nucleon cross section of 
both spin and scalar interactions as function of the WIMP mass are depicted below as scatter plots 
(Figs. 0-1). 

The use of the updated parameters (|^)-(^) does not change significantly the general distribution 
of points over the scatter plots as compared with calculations with earlier nucleon parameters 
f7%M used in iMITgl,^. 
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FIG. 1. Cross sections of spin-dependent and spin-independent interactions of WIMPs with proton 
and neutron. Filled triangles (light circles) correspond to relic neutralino density 0.1 < flxho < 0.3 
(0.025 < n^hl < 1). 

Scatter plots with individual cross sections of spin-dependent and spin-independent interactions 
of WIMPs with proton and neutron are given in fig. |^ as functions of the LSP mass. In the figure 
light circles correspond to cross sections calculated under the old assumption that 0.025 < ^x^o < 
1. Filled triangles give the same cross section but the constraint on the flat and accelerating 
universe is imposed by 0.1 < ^x^o ^ One can see that the reduction of the allowed domain 
for the relic density does not significantly affect spin-dependent and the spin-independent WIMP- 
nuclcon cross sections, i.e. restriction to a flat and accelerating universe weakly affects these cross 
sections. 

The different behavior of these cross sections with mass of the LSP can be seen from the plots. 
There is a more stringent lower bound for the spin-dependent cross section. It is at a level of 
10"'' pb, which is about an order of magnitude larger then the one presented in ]6^ , where for 
small tan {3 (tan /3 = 3, /i < 0, and small m^) the effect of a cancellation induced by the difference 
in signs between A„ and A^^s (@) was reported. Aside from the cancellation, the spin-dependent 
cross section peaks at about 10"'* pb and drops rapidly as increases down to 10"''"^® pb at 
nix ~ 600 GeV Q . We have checked that special consideration of the low tan f3 regime supplies 
us also with smaller cross section values for spin-dependent interactions, which do not enter in 
contradictions with . 

Such a cancellation was found also in scalar cross sections for tan/5 = 10 and /i < 0. In 
this case Higgs exchange is dominant. The cancellation in the mSUGRA is due to the cancellation 
between the up-type contribution (which is negative) and the down-type contribution, which is 
initially positive but decreasing, eventually becoming negative as we increase 
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FIG. 2. WIMP-nucleon cross section limits in pb for scalar interactions as function of the WIMP mass 
in GeV. Filled circles present our calculations with updated nucleon structure Q in "non-accelerating 
universe" with 0.025 < ^x^o < 1- Filled triangles give the same cross section but when 0.1 < fi^^/ig < 0-3. 
The contours (shaded area enclosed with solid curve) for allowed scalar WIMP-proton cross sections from 
§11 are also given together with some current (DAMA ||, HEIDELBERG-MOSCOwH, CDMS Q 
and HDMS prototype [pl|) and future experimental exclusion curves (HDMS projection [p^GENIUS-TF 
H, GENIUS 11 and CDMS H). 

If fig. H filled circles present our calculations when constraints due to an accelerating universe are 



not applied and as in 
but using as i,^,Pj67 



we hold 0.025 < fl^h^ < 1 



6? 



^1 



Filled triangles give the same cross section, 
the boundary 0.1 < il^h^ < 0.3. The contours for allowed 



scalar WIMP-proton cross sections from ||^,|6 8| a re also given together with some current (DAMA 



CDMS Q, HEIDELBERG-MOSCOW 11^, HDMS prototype |81 
perimental exclusion curves (HDMS [0, GENIUS-TF pf, GENIUS 



j) and future-expected ex- 
83i and CDMS [HI). This 



figure allows one to see the influence of the flat and accelerating universe on the distribution of 
WIMP-proton scalar cross section. The reduction left only 25% of points but nevertheless the dis- 
tribution of the remaining points differs only slightly from the one obtained with 0.025 < fi^/ig < 1. 
The models with very small cross sections as well as models with very large cross sections (in fact 
experimentally excluded) still persist. 
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FIG. 3. The same as in fig. but filled circles give cross sections for 0.1 < ^Ix^o < 0-3 tan/3 > 20. 
Filled triangle s g ive the same cross section, but when tan/3 > 40. The contour from (obtained for 
tan p < 10) ||,||| also given. 

One also can conclude that there is no contradiction between the resuh of obtained in the 

mSUGRA with minimal number of free parameters and our phenomenological scan, which also 
allows models with very small cross sections. 

While we have 1 < tan /3 < 50, the contour from [^6^ was obtained under the assumption 
that tan/3 < 10 to avoid some uncertainties in the treatment of radiative corrections in the 
renormalization-group evolution of the MSSM parameters which affect the relic density calcu- 
lations j6^. As noticed by many groups |58y24| ]6|, |35| , ^ , ^j22| , the scalar cross section of elastic 
WIMP-nucleon scattering increases with tan/3. As can be seen from fig. 1 of |^ tan/3 seems to be 
the only SUSY parameter with which the lower bound of the direct detection rate has the tendency 
to increase. The majority of the points at the scatter plots in fig. ^ are shifted to the domain of 
larger cross section with increase of tan/3. 

In general the increase of tan /3 effectively relaxes the /i constraint in mSUGRA (it allows /i to be 
smaller) and results in a non-negligible Higgsino component followed by significantly larger scalar 
cross section (see for example p|). 

We also respect as before [p8,E9Uq] non-universality of soft supersymmetry-breaking masses in 
the scalar and the gaugino sectors (see list of free parameters in (^)), resulting in larger cross 
sections, as noted in [ ]67| , p8[ . 

The spin-dependent and spin-independent WIMP-proton cross sections as functions of input 
parameters /i, mg. Ma, tan/3 are depicted in figs. ^ and ^. 
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FIG. 4. Cross sections of of WIMP-proton spin-dependent and spin-independent interactions as function 
of input parameters fi (upper panel) and tjiq (lower panel) obtained with 0.1 < fix^o < 0.3. 

There is no noticeable dependence of these scatter plots on the other free parameters from (|^), 
for which we therefore do not show scatter plots. One can see from fig. ^ the similarity of the 
scatter plots for spin-dependent and and scalar cross sections as functions of /j, and uiq . Decrease 
of both lower bounds of the cross sections with thq occur due to increase of masses of squarks, 
which enter the s-channel intermediate states. The only visible difference concerns more sharp 
lower bounds for the spin-dependent cross section. Both spin-dependent and spin-independent 
cross sections increase when |^| decreases, in agreement with [^,^||j3|]. It is not easy to trace 
the tendency in mSUGRA models because the parameter /i there is strongly constrained by the 
electroweak symmetry breaking condition (see, for example, [|68|). 

The increase of the scalar cross sections generally is connected with an increase of the Higgsino 
admixture of the LSP and increase of Higgsino-gaugino interference which enters this cross section 
|3^,^,^ . The reason of the Higgsino growth can be non- universality of scalar soft masses , 
variation of intermediate unification scale ]3^ ], or new focus point regime of supersymmetry pl| . 
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FIG. 5. Cross sections of WIMP-proton spin-dependent and spin-independent interactions as function 
of input parameters A/a (upper panel) and tan/3 (lower panel) obtained with 0.1 < fix'^o < 0.3. 

For example, as given in [Q, the smaller the intermediate scale Mj is, the larger the Higgsino 
components become. In particular, for Mj = 10^^ GeV the LSP is mainly Bino, the Higgs- 
neutralino-neutralino couplings are suppressed and therefore the cross sections are small. However 
for Mi ~ 10^^ GeV the Higgsino contributions become important and even dominant with the 
consequence of larger cross sections. It is also worth noting that, for any fixed value of Mj, the 
larger tan (3 is, the larger the Higgsino contributions become. 

Also it is claimed [ pT| that in the specific context of minimal supergravity (focus point regime), 
a cosmologically stable mixed gaugino-Higgsino state emerges as an excellent, robust dark matter 
candidate. The claim relies on recent arguments, that all squark and slepton masses can be taken 
well above 1 TeV with no loss of naturalness on the basis of a seemingly reasonable objective 
definition of naturalness . The mathematical basis of this result is the existence of focus points 
in renormalization group trajectories, which render the weak scale (i.e., the Higgs potential) largely 
insensitive to variations in unknown supcrsymmetry parameters^. While in these models the squark 



^The background of the approach can be questioned due to a possibility to shift the focus point for 
the Higgs mass parameter right to the GUT scale by means ofcappropriate choice of the initial condition for 
the top Yukawa coupling |8q] . Anyway from a phenomenological point of view the approach is interesting. 
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and slepton masses are unusually large, the electroweak gaugino and Higgsino particle masses are 
generically well below 1 TeV. The increase of the common soft scalar mass mo far beyond 1 TeV 
associated with decrease of below the gaugino masses, Mi,M2, leads to significant mixing 
between Higgsino and gaugino states accompanied by Higgs boson diagrams enhancement. A net 
result are again large scalar cross sections. 

The arguments presented above do not work in SUGRA . The LSP as Higgsino- like is almost 
excluded by LEP constraints [ p7| even if the assumptions of universal soft supersymmetry breaking 
are relaxed, and Higgsino dark matter is certainly excluded if universality is assumed, as is the case 
here. In addition to the LEP constraints, this is because the value of fj, is predicted as a function 
of mi/2 and mg, placing the LSP firmly in the Bino-like region. The same considerations exclude 
an LSP with mixed Higgsino /gaugino content. 

In the SUGRA framework of the elastic scattering cross sections, which are predicted for 
the LSP mass m-^ lie in a comparatively narrow band. This is essentially because the LSP is 
always mainly Bino-like, so its couplings do not depend greatly on other MSSM parameters such 
as mo. The principal causes of broadening are the uncertainties in the hadronic inputs and the 
possibilities of cancellations that may reduce the cross sections for some specific values of the 
constrained MSSM parameters . 

Figure H shows that while the spin cross section displays almost full insensivity to /i and Ma 
(Higgs bosons do not contribute) the scalar cross section possesses remarkable dependence on these 
parameters. The cross section rather quickly drops with growth of the CP-odd Higgs mass 
and increases with tan (3 in accordance with results of p8| , p4| J^ j35| , |3^ , ^ , p2| . 

The different tan (3- and Myi-dependence of spin-dependent and spin-independent cross section 
as well as general about-4-order-of-magnitude excess of spin-dependent cross section over spin- 
independent cross sections may be important for observations p7|,p8[. 



C. Role of the spin 

To be more definite with the statement claimed above, in fig. ^ we present a comparison of total 
spin-dependent versus total spin- independent event rates in ^'^Ge (spin= 9/2) — as representative 
and one of the most promising isotopes for future construction of high-sensitivity detectors. 

Figure || shows the weak dependence of the ratio on mass of the LSP with the mean value being 
approximately 0.01-0.1. There are very large and very small values for the ratio practically for 
any given mass of the LSP. The spin-independent (scalar) contribution obviously dominates in the 
domain of large expected rates in the Germanium detector (i? > 0.1 events/day/kg) as was obtained 
before (see, for example [Q). But as soon as the total rate drops down to i? < 0.01 evcnts/day/kg 
or, equivalently, the scalar neutralino-proton cross section becomes smaller than 10^^ ^ lO^^*' pb, 
the spin-dependent interaction may produce a rather non-negligible contribution to the total event 
rate. Moreover if the scalar cross section decreases further {a < 10~^^ pb), it becomes obvious that 
the spin contribution alone saturates the total rate and protects it (see lower bounds in fig. |^ and 
H) from decreasing below R « 10~^ 10"'' events/day/kg ||]. 

This observation could be quite important for experiments actually looking for direct detection of 
dark matter, but not only for exclusion plots. Indeed, while scalar cross sections governed mostly 
by Higgs exchange can be rather small (when Higgs masses remain large enough, for example in the 
Next-to-Minimal Supersymmetric Standard Model |Q ) the spin cross section can not be arbitrary 
small, because the mass of the Z boson, which gives the dominant contribution, is well defined, 
provided one ignores any possible fine-tuning cancellations [pSf . 

Therefore, if an experiment with sensitivity 10~^ ~ lO^^vents/day/kg fails to detect a dark 
matter signal, an experiment with higher sensitivity (and non-zero spin target) will be able to 
detect dark matter particles only due to the spin neutralino-quark interaction. 
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FIG. 6. Ratio of spin-dependent event rate to the spin-independent event rate in ^^Ge isotope as function 
of LSP mass (upper left), total (spin-dependent plus spin-independent) event rate (lower left) and scalar 
cross section of neutralino-proton interaction (lower right) obtained with 0.1 < 0,^h^ < 0.3. The vertical 
line gives the expected sensitivity of GENIUS |^^. In the region above the horizontal line the spin 
contribution dominates. The total event rate versus gaugino fraction of LSP P also given (upper right), 



IV. CONCLUSION 

Recent measurements in modern cosmology have changed the expected fraction of the cold dark 
matter in the universe, new results for the nucleon structure were obtained, a new generation of 
high-sensitivity experimental detectors are under consideration. All these changes stimulated us 
to recalculate our previous analysis concerning detection of cold dark matter. 

To this end we explored the MSSM parameter space at the weak scale where new accelerator 
and cosmological constraints are respected. We restrict the relic neutralino density to be in the 
rang e 0.1 < n^K^ < 0.3. We considered the variation of the spin-independent and spin-dependent 
WIMP-nucleon cross sections and of the expected event rate in ^"^Ge, with parameters of the 
MSSM, uncertainties of the nucleon structure and other theoretical assumptions like universality 
and coannihilation. 

The main results of the exploration can be summarized as follows. 

1 . The results of our updated calculations fall in general agreement with calculations performed 
in mSUGRA as well as with other less restrictive approaches, which allowed larger variation of the 
cross sections and detection rates. 

2. The use of the updated parameters of the nucleon structure does not change significantly 
the general distribution of points over the scatter plots as compared with calculations with earlier 
nucleon parameters. 
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3. In accordance with previous estimations we found that the coannihilation does not signif- 
icantly change our main conclusions. We understand that our estimation of the coannihilation 
effect is somewhat indirect, but in the effMSSM approach there is no stringent correlation between 
parameters, which sometimes makes the coannihilation channels inevitable. 

4. The new cosmological constraint on the relic neutralino density (due to flat and accelerating 
universe) which is numerically used in the form 0.1 < f^^^o < ^■'^ o^'^ approach does not 
significantly affect the resulting scatter plots for neutralino-nucleon elastic cross sections. 

5. To single out (mostly in the mSUGRA) theoretically a dominant contribution to the cross 
section or event rate one usually relies on the knowledge of the LSP composition. For example, 
as discussed through this paper, if the Bino fraction is large, then the cross section is small. 
Numerically the situation is less transparent. As seen from fig. ^ (upper right) the overwhelming 
majority of points (the region of highest point density) has P ^ 1, or very small Higgsino admixture 
and one should expect negligible event rate. Nevertheless this is not the case. There are a lot of 
points with sizable event rate for P 1. Therefore qualitative estimations of the dominance of 
the given contribution on the basis, for example, of large gaugino fraction of the LSP (P > 0.9) 
can be quantitatively not always correct. 

6. We notice that for targets with spin-non-zero nuclei it might be the spin- dependent interaction 
that determines the lower bound for the direct detection rate when the cross section of the scalar 
interaction drops below about 10^^^ pb. If this occurs the spin nuclear detectors would have notable 
advantage comparing with spinless detectors, or may become the only way to observe SUSY via 
direct dark matter detection. 

Finally we would like to stress again the fact, clearly seen from figs. |^ and that to reliably 
investigate the SUSY parameter space and therefore to have a chance to beat accelerator experi- 
ments in searching for (or discovery of) the new physics (supersymmetry) one needs a GENIUS-like 
detector. 
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